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TITLE OF THE INVENTION 
ELECTRIC NETWORK SIMULATING METHOD, SIMULATING 
APPARATUS, AND MEDIUM FOR STORING SIMULATION PROGRAM 
BACKGROUND OF THE INVENTION 

5 The present invention relates to a simulating 

method and apparatus for grasping the operating status 
of an electric network in advance in designing the 
electric network, and a computer-readable storage 
medium for storing a simulation program. 

10 As a conventional method of solving an electric 

network, for example, as shown in FIG. 15, to obtain 
voltages at nodes node^ to node3 connected to the 
respective elements in a circuit having resistive 
elements R^ and R2 connected in series with a current 

15 source I]_, the current and voltage formulas of this 

circuit are represented by matrix (1) as follows: 
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where G is the conductance as the reciprocal of 
20 a resistance R. 

When the above matrix is simply expressed as 

I = GV, it can be replaced with a mathematical problem 

for calculating a V matrix by giving I and G matrices. 

In the above matrix (1), G is LU-decomposed to rewrite 
25 matrix (1) into I = LUV. The LU matrix is comprised of 

lower and upper triangular matrices having elements 

represented by matrix (2): 



2 



Gil G12 G13 
G 2 i G 2 2 G 2 3 
G31 G 32 G33 



L n 0 0 
L 21 L 22 0 
L31 L 32 L 33. 



1 U12 U13 
0 1 U 2 3 
0 0 1 



• • (2) 



If a Z matrix satisfying I = LZ is present, 
LZ = LUV holds. z = UV is obtained by eliminating 
L from the right- and left-hand sides. If Z is known, 
V can be obtained. 

The general term of the L and U matrices can be 
expressed by equation (3), and the Z matrix can be 
expressed by equation (4). The V matrix as the final 
solution can be expressed by equation (5). A correct 
solution can be obtained by sequentially solving these 
equations . 

m-1 

.m = G im - X L i j u ji 

3 = 1 
m-1 

= ( G mj - 2 L mk u kj> / L mm 

k = l Jj = m + l,m + 2 / -' 



<!l - S L lk z k> / L ii 
k = l 



i = l,2,- 



Z u n-i,k v k 
k=n-i+l 



• (5> 



V n -i = Zn _ — 

0,1,2,-v n-1 

FIG. 16 is a flow chart of the prior art. 
Basically, the above contents are formed into an 
algorithm. It should be noted that the U matrix is 
calculated by a division using the diagonal elements of 
the matrix L. The diagonal elements of the matrix L 
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are L^, L22 / L 33r "'• Values (including zero) set in 
these diagonal elements are determined depending on the 
G matrix that describes the circuit. 

If the circuit is comprised of only resistive 
5 elements, and the resistances do not change, the value 

of the G matrix is also fixed. It is, therefore, 
possible to make the diagonal elements of the matrix L 
non-zero. If the diagonal elements become zero, the 
condition branches to change the order of the elements 

10 of the G matrix in FIG. 16. 

Assume that a resistive element has a variable R 
value. When the impedance value of this resistive 
element changes depending on the status of the circuit, 
that the diagonal elements of the L matrix are made 

15 non-zero is not guaranteed beforehand. In an actual 

algorithm, that the diagonal elements are made non-zero 
must always be monitored. 

If a value except non-zero but infinitely close to 
zero is set in the diagonal element, a division error 

20 occurs in a numerical calculation expressed by a finite 

number of bits, and a large calculation error may occur 
in the subsequent calculation results. In particular, 
when the differential term of the terminal voltage or 
current of an element is given as a denominator, for 

25 example, when a calculation item in which a value 

obtained by subtracting a voltage at time immediately 
preceding given time from a voltage at the given time 



becomes a denominator, a voltage applied across the 
element is not known beforehand, but known after a 
simulation. The value of this denominator is unknown. 
When a differential coefficient is zero, i.e., when the 
5 voltage value at the time immediately preceding the 

given time and the voltage value at the given time are 
zero or almost zero, the calculation result becomes 
very poor due to numerical errors. The subsequent 
calculations are meaningless in practice. An algorithm 

10 for always monitoring the number of effective digits 

even for non-zero values is required. 

The above problems are inevitably posed when a 
solution is obtained by solving a matrix. Even if the 
concept of the basic solution is simple, an algorithm 

15 for solving the above problems is required, resulting 

in redundant computation as a whole. 

BRIEF SUMMARY OF THE INVENTION 
It is object of the present invention to provide 
a simulating method and apparatus capable of reliably 

2 0 simulating a circuit operation with a relatively simple 

algorithm by regarding the circuit operation as 
movement of particles in a pipe, and a computer- 
readable storage medium which stores a simulation 
program. 

2 5 The present invention is an electric network 

simulating method comprising the steps of: after 
setting element cells representing electric functions 



of a plurality of circuit elements, intersection cells 
representing functions of electric wiring intersec- 
tions, and connection pipes representing connections 
between the element cells and the intersection cells, 
defining a current of an electric network as the number 
of particles moving in the connection pipe per unit 
time, and defining a voltage of the electric network as 
the number of particles present in the connection pipe; 
on the basis of the definitions in the defining step, 
setting beforehand, in units of element cells, a rule 
expressing an electric function of each of the circuit 
elements in accordance with a state of the connection 
pipe connected to the element cell, and setting 
beforehand, in units of intersection cells, a rule so 
that the numbers of particles present in the connection 
pipes connected to the intersection cell are equal 
to each other and a sum of the numbers of particles 
transferred at the intersection cell becomes zero; 
transferring particles between the element cell and the 
connection pipe and between the intersection cell and 
the connection pipe on the basis of the rules set in 
the setting step; and simulating the state of the 
electric network by updating the number of particles 
passing through a given connection pipe per unit time 
and the number of particles present in the given 
connection pipe and performing transfer and updating 
processes at least once. 



As described above, according to the simulating 
method of the present invention, the electric network 
is defined as the flow of particles, unlike the 
conventional method of simulating an electric network 
by solving simultaneous equations. For this purpose, a 
circuit is replaced with and redefined by element cells 
representing the electric functions of the circuit 
elements, the intersection cells representing the 
functions of the electric wiring intersections, and the 
connection pipes representing the connections between 
the element cells and the intersection cells. The rule 
representing the electric function of each circuit 
element is set in units of element cells in accordance 
with the state of the connection pipe. The particles 
are transferred between the element cells and the 
connection pipes and between the intersection cells 
and the connection pipe in accordance with this rule. 
The number and moving amount of particles in each 
connection pipe, i.e., the voltage and current can be 
solely determined. In principle, as described above, 
all the currents and voltages are redefined as the 
amounts and movements of particles. In addition, 
active elements such as a switch and semiconductor 
element can be equivalently replaced with basic 
elements such as a resistive element, voltage source, 
and current source. According to this idea, the 
relative amount of particles and the moving amount per 



unit time in a connected connection pipe are simulated. 
A quick accurate simulation result can be obtained 
for .a large electric network unlike the conventional 
simultaneous equation scheme in which errors occur in 
5 principle and a very long processing time is required. 

The present invention is not limited to the above 
simulating method, but can extend to a simulating 
apparatus for practicing the simulating method, and 
a storage medium storing a computer program which is 

10 installed in a computer apparatus to perform an 

equivalent simulating method. The simulating method 
and apparatus and the storage medium can achieve a high 
processing speed and high-precision simulation result, 
which cannot be obtained in the conventional method of 

15 solving simultaneous equations, according to the gists 

described above. 

Additional objects and advantages of the invention 
will be set forth in the description which follows, and 
in part will be obvious from the description, or may be 

2 0 learned by practice of the invention. The objects and 

advantages of the invention may be realized and 
obtained by means of the instrumentalities and 
combinations particularly pointed out hereinafter. 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING 

2 5 The accompanying drawings, which are incorporated 

in and constitute a part of the specification, 
illustrate presently preferred embodiments of the 
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invention, and together with the general description 
given above and the detailed description of the 
preferred embodiments given below, serve to explain the 
principles of the invention. 
5 FIG. 1 is a circuit diagram showing an electric 

circuit example according to an embodiment of the 
present invention; 

FIG. 2 is a view showing a model in which 
particles pass through pipes in the circuit example of 
10 FIG. 1; 

FIGS. 3A, 3B, and 3C are views showing an 
intersection of an electric circuit, and a model 
for referring to a connection pipe and a model for 
an intersection cell according to this embodiment; 
15 FIGS. 4A and 4B are views showing a model for 

a diode according to this embodiment; 

FIGS . 5A and 5B are circuit diagrams showing 
a capacitor and an equivalent model for the capacitor 
according to this embodiment; 
20 FIGS. 6A, 6B, and 6C are circuit diagrams showing 

an inductor, and equivalent models for the inductor 
according to this embodiment; 

FIG. 7 is a circuit diagram of a leakage 
transformer; 

25 FIG. 8 is a circuit diagram showing the equivalent 

model for the leakage transformer according to this 
embodiment; 



FIG. 9 is a circuit diagram showing a half -wave 
voltage resonance inverter; 

FIG. 10 is a circuit diagram showing the 
equivalent model for the half-wave voltage resonance 
inverter according to this embodiment; 

FIG. 11 is a flow chart showing a transient 
operation analysis process according to this 
embodiment ; 

FIG. 12 is a flow chart showing a periodic 
operation analysis process according to this 
embodiment ; 

FIG. 13 is a block diagram showing the arrangement 
of a simulating apparatus according to this embodiment; 

FIG. 14 is a view showing an output example in the 
simulating apparatus; 

FIG. 15 is a circuit diagram for explaining 
a conventional method of solving an electric network; 

FIG. 16 is a flow chart showing the arithmetic 
process of the electric network; and 

FIG. 17 is a table showing a cell list according 
to this embodiment. 

DETAILED DESCRIPTION OF THE INVENTION 

A preferred embodiment of the present invention 
will be described below with reference to the 
accompanying drawings . 

FIG. 1 shows a closed loop circuit in which 
a resistor R is connected to a current source 1. 
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Voltages VI and V2 appear at two terminals of the 
resistor R. The current source 1 supplies a current i. 
A particle model equivalent to this circuit is shown in 
FIG. 2. 

5 In this particle model, the current source 1 is 

replaced with a pump P. This pump P receives predeter- 
mined particles per unit time from a lower pipe 2 and 
outputs the particles to an upper pipe 3. A function 
corresponding to the resistor R is expressed as a 

10 restrictor for particles passing from the upper pipe 3 

to the lower pipe 2. A unit volume is used for simply 
describing the internal volumes of the pipes 2 and 3 . 
The number of particles in each of the pipes 2 and 3 
directly represents the particle density. 

15 When this particle density is regarded as a 

voltage, and the number of particles moving per unit 
time is regarded as a current, the restriction of 
the particles at the restrictor corresponding to the 
resistor R can be expressed as a physical amount. 

20 More specifically, for example, when the current source 

1 has 1 A, and the resistor R has 100 Q, a potential 
difference across the resistor R is 100V. In this 
case, 1 A is expressed by 1,000 particles, and the 
degree of restriction of the resistor R is directly 

25 expressed as 100. 

In a balanced state, the difference corresponding 
to 100,000 (= 1,000 X 100) particles is produced across 



the resistor R. The number of particles is not 
necessarily positive, but may be zero or negative. 
The final necessary result is the difference in the 
number of particles across the resistor R. 

FIG. 3A shows wiring connections in an electric 
circuit. An intersection in a general electric circuit 
is a simple connection and is not defined as an 
element. All nodes at the intersections are defined 
as a and are not discriminated from each other in the 
circuit. 

To the contrary, as shown in FIG. 3C, an intersec- 
tion in a particle model is regarded to have the same 
function as that of an element, and all nodes connected 
to an intersection cell are different from each other. 
That is, three nodes a, b, and c are present. The 
number of particles present in three pipes connected to 
this intersection cell are averaged, and the sum of the 
numbers of particles moving through the pipes at the 
time of averaging is set to zero. 

If only a voltage and current are referred to 
at an intersection, the voltage and current can be 
expressed as a combination of a given pipe through 
which particles are moving and a pipe for referring to 
the given pipe. 

FIGS. 4A and 4B show an element cell imitating 
a diode function. This element cell has two terminals 
having polarities. For example, as shown in FIG. 3A, 



assume that the cathode and anode are defined to have 
positive and negative polarities , respectively. If the 
number of particles in the cathode-side pipe is larger 
than that in the anode-side pipe, no particle movement 
occurs. However, as shown in FIG. 3B, when the cathode 
and anode are negative and positive, respectively, the 
particles smoothly move from the anode to the cathode 
with an increase in the number of particles. 

A description of an element cell having this 
function allows obtaining diode characteristics as 
the rectification function. If more strict diode 
characteristics are required, i.e., if a junction 
capacitance, junction potential, temperature coeffi- 
cient, and breakdown voltage are also included, a cell 
may be defined as an arbitrary function including these 
factors. If it is redundant to describe the diode as 
one element cell, a diode may be described in detail as 
a combination of basic linear elements. 

A capacitor element will be modeled. A capacitor 
can be expressed as a current integration formula of a 
terminal voltage V, as indicated by equation (6) below. 
If this is regarded as the sum of regions segmented 
with finite shortest time intervals, the integration 
symbol can be expressed as a £ function. This can be 
expressed as the primary term between the current and 
a terminal voltage V old of the immediately preceding 
time . 
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... (6) 



The above equation as a linear equation between 
the current and voltage has a voltage source as V Q id 
5 and a resistive component as dt/C. When this 

arrangement is decomposed into basic linear elements 
and equivalently expressed, the capacitor C shown in 
FIG. 5A- can be expressed as a series connection of 
a resistive element dt/C and a voltage source V old , as 
10 shown in FIG. 5B. 

A calculation for an arrangement including this 
capacitor is performed at given time t, and its 
operating state at time immediately preceding the given 
time is already confirmed. In the calculation at given 
15 time t, V old can be regarded as a voltage source having 

a fixed value. The shortest time section dt of the 
resistive element is known, and the capacitance of the 
capacitor C is also known. These values are also fixed 
values. Therefore, the behavior of the capacitor at 
2 0 time t can be expressed by these basic elements. 

An inductor element will be modeled below. As 
represented by equation (7), the left-hand side may be 
described with v as in the capacitor. in this case, 
the right-hand side becomes a current differential 
25 term. Equation (7) can be solved, but a current value 

in an insufficient convergent state is inaccurate. 
It may often be unpreferable to apply the differential 
to an inaccurate current. 



As indicated by equation (8), the current is 
defined on the left-hand side, and the voltage 
integration value is defined on the right-hand side. 
The current at time t is the sum of the current value 
L/Dt X ioXcj flowing at time immediately preceding time 
t and the current flowing through the resistor, as 
shown in FIG. 6B. If this is expressed using basic 
elements, the inductor L shown in FIG. 6A is expressed 
as a parallel connection between the resistive element 
L/dt and the current source i Q ld' as shown in FIG. 6C. 



di L . L . 
dt dt dt ° 



ifvdt 

L J 



^2> = iold + ^v ... (8) 



FIG. 7 shows the general notation of a leakage 
transformer T. In this transformer T, two inductors, 
i.e., Lp and Ls are magnetically coupled. The degree 
of coupling is expressed by a coupling coefficient k. 
The equations holding on the primary and secondary 
sides of this transformer T are given by equations (9), 
(10), and (11) 

V p = L p ^P + M ps ^ -(9) 

V s = L s + M ps •••(10) 

S S dt PS dt 

Mps = k^L p L s ••• (11) 

where Mp S is the transinductance . 

The modification of equations (9), (10), and (11) 
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in accordance with the same idea as in the inductor 
yields equations (12) and (13): 

- ■ dt ^-V S -(12) 



ippast + Lp(1 _ k2) V P - LpLs _ M 2 s 
dt Mdt 



ls ^ St L S (1 - k2, V * - LpLg _ Mps 



Vp •••(13) 



where i ppa st and ispast are currents flowing through 
the primary and secondary sides at time immediately 
preceding time t, and V p and V s are terminal potential 
differences on the primary and secondary sides at 
time t. 

These equations are implemented with element cells 
as a model shown in FIG. 8. That is, this model can be 
expressed as current source cells Ip pa st and I Spast^ 
resistor cells R-tpI and R TS1' R TM1' and R TM1' and 
intersection cells. 

How to recognize this circuit diagram will be 
described below. 

FIG. 9 shows a one-resonator type half -wave 
voltage resonance inverter, in this inverter, a 
capacitor C2 is connected to a DC power supply VI via 
a parallel circuit of a diode Dl and capacitor CI. 
A series circuit of a switch SI and a parallel circuit 
of a capacitor C3 and a primary side LI of the leakage 
transformer T is connected to the capacitor C3 . 
A diode D2 is connected in parallel to the switch SI 
in an opposite polarity, thereby forming a half -wave 
switch circuit. 
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A load circuit Rl is connected to a secondary side 
L2 of the transformer T via a capacitor C4. A series 
circuit of a resistor R2 and capacitor C5 is connected 
in parallel to the load circuit Rl . 

5 This inverter is powered by the DC power supply 

VI. The half -wave switch circuit made up of the switch 
SI and diode D2 turns on/off a parallel resonator made 
up of the capacitor C3 and the primary side LI to 
continue parallel resonance. Part of the resonance 

10 energy is transmitted to the secondary side using the 

leakage component of the transformer T to drive the 
load circuit Rl. 

This inverter can be replaced with the element 
cells described above as an equivalent cell circuit 

15 diagram shown in FIG. 10. More specifically, the 

capacitor C2 is connected in series with a resistive 
element R C 2 and voltage source V C 2 . The capacitor C3 
is connected in series with a resistive element R C 3 and 
voltage source V C 3 . The capacitor C4 is connected in 

2 0 series with a resistive element Rc4 and voltage source 

Vc4- The capacitor C5 is connected in series with 
a resistive element R C 5 and voltage source V C 5 . The 
leakage transformer T is expressed as a combination of 
current source cells ItpI' and I TS1' resistor cells 

2 5 RtpI' r TS1' r TM1' and R TM1' and intersection cells X4 , 

X5, X6, and X7 . 

Reference symbols Yl and Y2 in FIG. 10 denote 



intersection cells each at which three pipes intersect. 
Reference symbols XI, X2 , •*•, X8 denote intersection 
cells each at which four pipes intersect. The internal 
processes of these pipes are the same although the 
5 numbers of intersecting pipes are different. The 

switch SI is expressed as a turn-on/off function as 
in the diode. The ON /OFF state of the switch SI is 
determined by a voltage (number of particles) applied 
to the control terminal. A pulse source PI is 

10 connected to the control terminal, and its reference 

potential is given as a switch control voltage obtained 
by making a reference cell (i.e., a cell for referring 
to the number of particles in a pipe to be referred to) 
refer to the number of particles at a node 15 and 

15 adding the number of particles corresponding to a 

predetermined pulse voltage to the number of particles 
referred to by the reference cell. The node names are 
identical at the intersections in FIG. 9. However, 
since nodes are newly assigned in FIG. 10, nodes 0 to 

20 39 are set in FIG. 10, while nodes 0 to 7 are set in 

FIG . 9. 

An array shown in FIG. 17 is defined as informa- 
tion for inputting such an equivalent cell arrangement. 
NME$ in FIG. 17 is an array for storing the name of an 
25 element cell or intersection cell. The first character 

of the name represents the attribute of the cell. 
The index of the array starts from "1", and the number 
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of indices is the number of cells. Since the total 
number of cells is 32, indices are "1" to "32". 

Cells (element cells and intersection cells) 
required in the present invention are listed in 
5 FIG. 17. 

According to the cell list in FIG. 17, NP and 
NM represent arrays for storing node numbers of nodes 
to which a cell is connected. When the number of 
terminals of an element is three or more, arrays N3 , 

10 N4, Nn are used following the arrays NP and NM. 

When an element has polarities, a polarity description 
becomes important. The polarity has particularly no 
meaning for an element having no polarity, such as a 
resistor. However, for convenience, the same notation 

15 is employed. DTA represents an array for storing 

unique information of a cell. A voltage value is set 
in a voltage source, a resistance value is set in 
a resistive element, and no value is set in a cell 
such as an intersection cell that does not require 

20 particular unique information. When an electric 

element such as a capacitor, inductor, or leakage 
transformer, which is expressed as a combination of 
more basic equivalent cells, the respective equivalent 
cells describe unique information of the cell. In this 

25 case, the unique information is obtained from the 

above-described equivalent equations. 

A variable new_i described last loops the process 
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the number of cells. The number of particles present 
in a pipe is stored in an array V(i), and the number of 
particles passing per unit time is stored in an array 
II (i), where i is not the cell number but the pipe node 
5 number. The number of arrays V or II is equal to the 

number of pipes. If an initial value is required in 
analysis of an arrangement including an inductor and 
capacitor, the values obtained in equations (6) to (8) 
described above are set as initial data of the arrays V 
10 and II. 

A practical example of a program for transferring 
particles in cells will be described below. For 
example, a process for a resistor cell starting from R 
is described as follows. 
15 Case"R" 

v_temp=V(NP(i) ,t_cell)-V(NM(i) ,t_cell) 
i_temp=V_temp/DTA( i ) 

V(NP(i) ,t_cell)=V(NP(i) , t_cell ) -i_temp 
V(NM(i) ,t_cell)=V(NM(i) ,t_cell)+i_temp 

2 0 wherein V is the array for storing the number of 

particles. If an array is defined as V(i,t) including 
the time dimension, i represents the node number 
assigned to a pipe. The number i is not a number 
corresponding to the loop order of the variable new_i 

2 5 defined by the array NM$(i), and t represents the 

number of steps corresponding to calculation time. 
To obtain a steady solution, the t term is unnecessary. 
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When transient operation analysis is required, t is 
updated shortest time interval. 

V(NP(i) ,t_cell) represents the number of 
particles for the number of steps t_cell in a pipe 
5 represented by a positive node array NP(i). 

Similarly, V(NM(i) ,t_cell) is the number of particles 
in a pipe represented by a negative node array NM(i). 
A difference v_temp between the above numbers of 
particles is the difference between the numbers of 

10 particles in pipes at two terminals of the resistive 

cell. This difference represents the potential 
difference. Since the resistance value is input in the 
DTA array for a resistor cell, a particle flow corre- 
sponding to the current to be flowed can be expressed 

15 as i_temp=V_temp/DTA(i) . In the resistor cell, 

particles flow from the positive terminal, and the 
particles equal in number to the input particles are 
output from the negative terminal. The number of 
particles corresponding to i_temp is added to or 

2 0 subtracted from the numbers of particles in the pipes 

at the two terminals of the resistor cell. The 
foregoing is the cell function as that of the resistor. 

The diode cell is a cell having a rectification 
function of flowing a current in the forward direction 

25 but flowing no current in the opposite direction. 

The diode cell defines a voltage corresponding to the 
p-n junction potential when the current flows in the 
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forward direction. 
Case"D" 

V_temp=V ( NP ( i ) ,t_cell)-V(NM(i) ,t_cell) 
IF V_temp>0.7 Then 
5 V(NP(i) ,t_cell)=V(NP(i) ,t_cell)-V_temp/2 

V(NM(i) ,t_cell)=V(NM(i) / t_cell)+V_temp/2 
End if 

The potential difference V_temp applied to the 
diode cell in the same manner as described above is 

10 calculated. When NP(i) and NM ( i ) are defined as the 

anode and cathode sides of the diode, respectively, 
V_temp is the potential difference when viewed from the 
anode side to the cathode side and corresponds to the 
forward voltage. Assume that the particles move when 

15 the forward voltage is 0.7V or more, and the number of 

moving particles is given as V_temp/2. The number of 
particle at one terminal of the element must be equal 
to that at the other terminal of the element. A 
description is made such that the number of particles 

2 0 V_temp/2 is subtracted from that on the anode side and 

added to that on the cathode side. When the forward 
voltage is 0.7V or less, the diode is cut off, and no 
particles move through the diode. In this case, no 
description need be made for particle movement. 

25 An intersection cell having, e.g., three terminals 

can be described as follows. An intersection in 
an electric circuit must receive and output a common 
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voltage and have the total sum of input/output currents 
as zero. These two conditions are necessary in solving 
simultaneous equations. The function to be described 
for the intersection cell is to satisfy the condition 
5 that the numbers of particles between the connected 

pipes are equal to each other. When this intersection 
cell functions, the total sum of input/output particle 
flows (numbers of particles per unit time) consequently 
becomes zero in the balanced state. In the particle 
10 model, obviously, the number of particles (voltage) 

is the basic physical amount, while the particle flow 

(current) is a subordinate physical amount. 

Case"Y" 

V_temp= 

15 (V(NP(i) ) ,t_cell)+V(NM(i) ,t_cell)+V(N3(i) ,t_cell)/3 

V(NP(i) ,t_cell)=v_temp 
V ( NM ( i ) , t_cell ) =V_temp 
V(N3 (i) ,t_cell)=V_temp 

A current source cell sets a status in which a 
2 0 constant current flows in a steady state. The number 

of particles per unit time is input and output to the 
two end pipes . Since the number of times the scan loop 
of the cell is performed per unit time is proportional 
to the current value, a value set in the DTA array is 
25 a coefficient multiple of the current value. The most 

understandable method is to set the current value 
itself in the DTA array and perform the scan loop once 
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per unit time. 
Case" I" 

V(NP(i) ,t_cell)=V(NP(i) ,t_cell)+DTA(i) 
V(NM(i) ,t_cell)=V(NM(i) ,t_cell)=DTA(i) 
5 An analysis algorithm using these cells will now 

be described. 

FIG. 11 is a flow chart showing an example of 
transient operation analysis. When the program runs, 
an array declaration is performed in SI. That is, 
10 a memory area for necessary arrays and variables is 

assured. In S2 , information for a circuit to be 
analyzed is input. The input circuit information is 
circuit data of a general electric circuit. 

In S3, equivalent cell information is extracted. 
15 In S4, pipe information is extracted. More specifi- 

cally, since the circuit data of the electric circuit 
includes data that cannot be directly expressed as 
element cells, the circuit data are developed as 
circuit information of equivalent cells. In this case, 
20 intersections are replaced with intersection cells, and 

node numbers shared by the intersection positions are 
replaced with those of the intersection cells. The 
node numbers are defined as entirely different node 
numbers. These node numbers are used to identify 
25 pipes. At this moment, a cell list in FIG. 17 is 

almost finished. 

In S5, initial values are set. More specifically, 
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the initial electric state for executing a simulation 
is set. For example, the initial value for a capacitor 
is an initial voltage, and the initial value for an 
inductor is an initial current. The initial voltage 
5 indicates that the initial number of particles present 

in a pipe is determined, i.e., the initial value of the 
V array is determined. The initial current indicates 
that the number of particles flows into or out from an 
inductor per initial unit time is determined as the II 

10 array beforehand. 

A loop structure will be described. The "transfer 
of particles in each cell" in this loop structure 
indicates that particles are transferred in the cells 
shown in FIG. 17 in accordance with the attributes 

15 and pieces of unique information of these cells. 

This transfer is performed for all the cells. 

The first loop (S7) is to simulate the steady 
state of the circuit at given time. In the first loop, 
the states of the particles are not balanced as a whole 

20 at the beginning. When transfer (S8) is repeated 

several times, the numbers of particles and particle 
movements in the pipes converge. When a predetermined 
convergence condition is satisfied, the state at the 
given time is regarded to be confirmed, and the first 

25 loop is ended. The convergence is determined such that 

in the first loop, the differences between previous and 
current values of all or important parts of the V and 



II arrays are smaller than the convergence error (S9) 
and are compared with predetermined convergence 
condition . 

When the first loop is ended, the time advances 
a shortest time interval (S9), and the flow returns to 
the first loop again- This repetition is performed in 
the second loop (S6). This loop is repeated from time 
0 to desired time (Sll) and ended (S12). When the time 
advances the shortest time interval, the following 
operation is performed. For example, equivalent 
conversions shown in FIGS. 5A and 5B and FIGS. 6A and 
6B are performed for the capacitor and inductor. 
For example, in the inductor equivalent model shown in 
FIG. 6C, a new impedance obtained upon the lapse of the 
shortest time interval is calculated by regarding as 
each II array the total sum of the numbers of particles 
per unit time flowing in the equivalent model at time 
preceding the shortest time interval. This impedance 
is newly set as the unique information of the current 
source i 0 ld* Tne unique information of the voltage 
source of the equivalent model (FIG. 5B) for the 
capacitor is similarly calculated from the value of the 
V array at time preceding the shortest time interval. 

In this transient operation analysis, convergence 
in the first loop is important and must be strictly 
determined. If this determination is less strict, the 
result adversely affects the subsequent operations on 
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the time basis. 

FIG. 12 is a flow chart showing an example of 
periodic operation analysis. The basic flow (S21 
to S24, S27 to S32, and S35) is the same as in the 
5 transient operation analysis described above, except 

that the second loop in the transient operation 
analysis is determined on the basis of time, while the 
second loop is ended at the end of the period in the 
periodic operation analysis (S33), and the third loop 

10 is newly installed in the periodic operation analysis. 

In the periodic operation analysis, a one-period 
simulation is complete at the end of the second loop. 
This indicates that one-period transient operation 
analysis is performed for the given initial condition. 

15 In an actual circuit as well, the periodic operation 

stabilizes to obtain a steady state several cycles 
after the start with the initial condition. This also 
applies to the simulation. To obtain the steady state 
in the periodic operation, the third loop (S34 and S26) 

2 0 must be repeated several times. 

The above operation is equivalent to repetition 
of transient operation analysis, but the periodic 
operation analysis has an advantage in that the 
convergence condition in the first loop can be less 

25 strict at the beginning. More specifically, the 

process for obtaining the steady solution of the period 
need not always trace a transiently correct process. 
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Rather, the third loop is repeated a larger number of 
times to make an element having a large time constant 
converge . 

FIG. 13 is a block diagram showing the arrangement 
5 of a simulating apparatus for implementing the above 

simulation. The simulating apparatus is comprised of a 
simulation unit 11, input device 12, and output device 
13. An example of the input device 12 is a keyboard 
for directly inputting circuit connection information 

10 with key operations, or a storage medium such as a 

floppy disk or magnetooptical disk for storing circuit 
connection information in the form of a file and 
loading the stored circuit connection information to 
the simulation unit 11. An example of the output 

15 device 13 is a display for displaying a simulation 

result on the monitor screen or a printer for printing 
it out. 

The simulation unit 11 is comprised of an input 
unit 14, arithmetic unit 16, and output unit 17. 

2 0 The input unit 14 converts the circuit connection 

information from the input device 12 into information 
suitable for internal processing and receiving the 
converted information. The arithmetic unit 16 performs 
the following arithmetic operations. That is, the 

25 arithmetic unit 16 writes the information from the 

input unit 14 in a circuit information storage area 151 
of a storage unit 15, creates an equivalent cell model 
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on the basis of the written information, and stores the 
model information in a cell array storage area 152. 
The arithmetic unit 16 also sequentially stores the 
number of particles in a pipe connected to each cell 
5 and the number of particles moving per unit time in a 

particle status storage area 153 and sets an initial 
value in the storage area 153. The output unit 17 
converts the arithmetic result of the arithmetic unit 
16 into a format suitable for the output device 13. 

10 An output example on the monitor screen upon 

simulating a circuit operation using the simulating 
apparatus having the above arrangement is shown in 
FIG. 14. This output example shows the simulation 
result obtained using the circuit arraignments shown in 

15 FIGS. 9, 10, and 17 in accordance with the flow chart 

of periodic operation analysis shown in FIG. 12. The 
constants of the circuit elements are listed and output 
on the left side of the screen, while the node voltage 
and current waveforms that are desired at the switching 

2 0 period are displayed on the right side of the screen. 

The conditions for convergence determination are 
displayed in the four windows on the upper right 
position of the screen. 

Using this simulating apparatus, the operating 

25 status of the circuit can be easily grasped, and 

circuits can be efficiently designed. More specifi- 
cally, since circuit operation is basically defined as 
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the movement of particles in each pipe, the arithmetic 
contents are additions or subtractions. Since no 
divisions are performed, an overflow as the essential 
factor of an algorithm does not occur. Calculations 
5 are not interrupted and simulations within 

predetermined time ranges can be performed regardless 
of combinations of the circuit arrangements and circuit 
elements . 

Note that although the additions and subtractions 
10 are performed for pipes, divisions may be included 

depending on the element functions of element cells. 
For example, in a description of an element cell 
representing a resistive element, a moving amount of 
particles corresponding to a current is calculated by a 
15 division using the resistance value. If the resistance 

value is set to zero, an overflow occurs. However, the 
resistance value is known beforehand, and appropriate 
processing can be simply performed. That is, when 
element cells are properly described, the arithmetic 
20 operations can be separated from the simulating method 

itself . 

The storage medium such as a floppy disk or 
magnetooptical disk can also store the simulating 
processing program itself executed by the simulation 
2 5 unit 11 in addition to the circuit connection informa- 

tion. The simulation unit 11 loads the simulation 
processing program from the storage medium and then 



loads the circuit connection information to perform 
circuit operation simulations. 

As has been described above , according to the 
present invention, circuit operations are defined as 
movements of particles in pipes. There can be provided 
an electric network simulating method capable of 
reliably simulating circuit operations with a 
relatively simple algorithm. 

In addition , according to the present invention, 
circuit operations are defined as movements of 
particles in pipes. There can be provided an electric 
network simulating apparatus capable of reliably 
simulating circuit operations with the relatively 
simple algorithm. 

According to the present invention, there can also 
be provided a storage medium that stores a simulation 
program capable of reliably simulating circuit 
operations with the relatively simple algorithm. 
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WHAT IS CLAIMED IS: 
1. An electric network simulating method 
comprising the steps of: 

defining element cells representing electric 
5 functions of a plurality of circuit elements and 

connection pipes representing wiring lines for 
connecting the circuit elements, defining an electric 
network current as the number of particles moving in 
the connection pipe per unit time, and defining an 
10 electric network voltage as the number of particles 

present in the connection pipe; 

on the basis of definitions in the defining step, 
setting beforehand, in units of element cells, a rule 
for expressing an electric function of each of the 
15 circuit elements in accordance with a state of the 

connection pipe connected to each of the element cells; 

transferring particles between the element cell 
and the connection pipe in accordance with the set 
rule; and 

20 simulating the state of the electric network by 

updating the number of particles passing through a 
given connection pipe per unit time in the transferring 
step and the number of particles present in the given 
connection pipe, and performing transfer and updating 

25 processes at least once. 

2. An electric network simulating method 
comprising the steps of: 
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after setting element cells representing electric 
functions of a plurality of circuit elements, 
intersection cells representing functions of electric 
wiring intersections, and connection pipes representing 
5 connections between the element cells and the 

intersection cells, defining a current of an electric 
network as the number of particles moving in the 
connection pipe per unit time, and defining a voltage 
of the electric network as the number of particles 

10 present in the connection pipe; 

on the basis of the definitions in the defining 
step, setting beforehand, in units of element cells, 
a rule expressing an electric function of each of the 
circuit elements in accordance with a state of the 

15 connection pipe connected to the element cell, and 

setting beforehand, in units of intersection cells, 
a rule so that the numbers of particles present in 
the connection pipes connected to the intersection cell 
are equal to each other and a sum of the numbers of 

2 0 particles transferred at the intersection cell becomes 

zero; 

transferring particles between the element cell 
and the connection pipe and between the intersection 
cell and the connection pipe on the basis of the rules 
25 set in the setting step; and 

simulating the state of the electric network by 
updating the number of particles passing through a 
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given connection pipe per unit time and the number of 
particles present in the given connection pipe in the 
transferring step and performing transfer and updating 
processes at least once. 
5 3 . An electric network simulating method 

according to one of claims 1 and 2, wherein the setting 
step includes the step of 

when a given one of the circuit elements is 
a current source, setting a rule for extracting the 

10 number of particles corresponding to a current value 

per unit time from one of two connection pipes 
connected to an element cell expressing the given 
circuit element and giving the number of particles 
equal in number to the number of extracted particles 

15 to the other one of the two connection pipes. 

4. An electric network simulating method 
according to one of claims 1 and 2, wherein the setting 
step includes the step of 

when a given one of the circuit elements is a 
2 0 voltage source, setting a rule for making a difference 

between the number of particles in one of two connec- 
tion pipes connected to an element cell expressing the 
given circuit element and the number of particles in 
the other one of the two connection pipes equal to the 
25 number of particles corresponding to a voltage of the 

voltage source . 

5. An electric network simulating method 
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according to one of claims 1 and 2, wherein the 
defining step includes the step of 

when a given one of the element cells has 
nonlinearity as a function of time, defining the given 
5 circuit element as a combination of an element cell 

for a resistive element and one of an element cell 
expressing a current source and an element cell 
expressing a voltage source, the combination expresses 
linearity equivalent to a behavior of the given circuit 

10 element at given time; and 

the setting step includes the steps of 
when a certain one of the circuit elements is 
a current source, setting a rule for extracting the 
number of particles corresponding to a current value 

15 per unit time from one of two connection pipes 

connected to an element cell expressing the certain 
circuit element and giving the number of particles 
equal in number to the number of extracted particles 
to the other one of the two connection pipes, and 

20 when a specific one of the circuit elements is a 

voltage source, setting a rule for making a difference 
between the number of particles in one of two connec- 
tion pipes connected to an element cell expressing the 
specific circuit element and the number of particles in 

2 5 the other one of the two connection pipes equal to the 

number of particles corresponding to a voltage of the 
voltage source. 
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6. An electric network simulating method 
according to one of claims 1 and 2, wherein the setting 
step includes the step of 

when a given one of the circuit elements has 
5 an impedance characteristic discontinuously changing, 

preparing a plurality of rules for the element cell for 
expressing the given circuit element and selecting one 
of the plurality of rules in accordance with the state 
of the connection pipe connected to the element cell. 

10 7. An electric network simulating method 

according to claim 5, wherein the transferring step 
and the simulating step include the step of 

simulating the state of each element cell at 
initial time so as to simulate a transient phenomenon 

15 of the given circuit element having nonlinearity as 

a function of time, simulating a behavior of the 
nonlinear element at an operating point advancing by 
a shortest time interval, by changing each parameter of 
a combination of the element cells having functions 

2 0 equivalent to the element cells, and simulating the 

transient phenomenon by repeating the change in 
parameter every time the shortest time interval has 
elapsed. 

8. An electric network simulating method 
25 according to claim 6, wherein the transferring step 

and the simulating step include the step of 

simulating a behavior of each element cell at 
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initial time so as to simulate a transient phenomenon 
of the given circuit element having the impedance 
characteristic discontinuously changing, simulating 
a behavior of the nonlinear element at an operating 
5 point advancing a shortest time interval by executing 

the transferring step in accordance with the rule 
selected in accordance with the state of the connection 
pipe connected to the element cell, and simulating the 
transient phenomenon by repeating the simulating steps 
10 every time the shortest time interval has elapsed. 

9. An electric network simulating apparatus 
comprising: 

means for defining element cells representing 
electric functions of a plurality of circuit elements 

15 and connection pipes representing wiring lines for 

connecting the circuit elements, defining an electric 
network current as the number of particles moving in 
the connection pipe per unit time, and defining an 
electric network voltage as the number of particles 

2 0 present in the connection pipe; 

means for setting beforehand, on the basis of 
definitions in the defining means, in units of element 
cells, a rule for expressing an electric function of 
each of the circuit elements in accordance with a state 

2 5 of the connection pipe connected to each of the element 

cells ; 

means for transferring particles between the 
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element cell and the connection pipe in accordance with 
the set rule; and 

means for simulating the state of the electric 
network by updating the number of particles passing 
5 through a given connection pipe per unit time in the 

transferring means and the number of particles present 
in the given connection pipe, and performing the 
transfer means and updating process at least once. 
10. An electric network simulating apparatus 

10 comprising: 

means for, after setting element cells 
representing electric functions of a plurality of 
circuit elements, intersection cells representing 
functions of electric wiring intersections, and 

15 connection pipes representing connections between the 

element cells and the intersection cells, defining 
a current of an electric network as the number of 
particles moving in the connection pipe per unit time, 
and defining a voltage of the electric network as the 

2 0 number of particles present in the connection pipe; 

means for setting beforehand, on the basis of the 
definitions in the defining means, in units of element 
cells, a rule expressing an electric function of each 
of the circuit elements in accordance with a state of 

2 5 the connection pipe connected to the element cell, and 

setting beforehand, in units of intersection cells, 
a rule so that the numbers of particles present in the 
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connection pipes connected to the intersection cell 
are equal to each other and a sum of the numbers of 
particles transferred at the intersection cell becomes 
zero; 

5 means for transferring particles between the 

element cell and the connection pipe and between the 
intersection cell and the connection pipe on the basis 
of the rules set in the setting means; and 

means for simulating the state of the electric 

10 network by updating the number of particles passing 

through a given connection pipe per unit time and the 
number of particles present in the given connection 
pipe in the transferring means and performing transfer 
and updating processes at least once. 

15 11. An electric network simulating apparatus 

according to one of claims 9 and 10, wherein the 
setting means includes 

means for setting, when a given one of the circuit 
elements is a current source, a rule for extracting the 

2 0 number of particles corresponding to a current value 

per unit time from one of two connection pipes 
connected to an element cell expressing the given 
circuit element and giving the number of particles 
equal in number to the number of extracted particles 

25 to the other one of the two connection pipes. 

12. An electric network simulating apparatus 
according to one of claims 9 and 10, wherein the 
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setting means includes 

means for setting, when a given one of the circuit 
elements is a voltage source, a rule for making a 
difference between the number of particles in one of 
5 two connection pipes connected to an element cell 

expressing the given circuit element and the number of 
particles in the other one of the two connection pipes 
equal to the number of particles corresponding to 
a voltage of the voltage source. 

10 13. An electric network simulating method 

according to one of claims 9 and 10, wherein the 
defining means includes 

means for defining, when a given one of the 
element cells has nonlinearity as a function of time, 

15 the given circuit element as a combination of an 

element cell for a resistive element and one of an 
element cell expressing a current source and an element 
cell expressing a voltage source, the combination 
expresses linearity equivalent to a behavior of the 

2 0 given circuit element at given time; and 

the setting means includes 

means for setting, when a certain one of the 
circuit elements is a current source, a rule for 
extracting the number of particles corresponding to a 
25 current value per unit time from one of two connection 

pipes connected to an element cell expressing the 
certain circuit element and giving the number of 
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particles equal in number to the number of extracted 
particles to the other one of the two connection pipes, 
and 

means for setting, when a specific one of the 
5 circuit elements is a voltage source, a rule for making 

a difference between the number of particles in one 
of two connection pipes connected to an element cell 
expressing the specific circuit element and the number 
of particles in the other one of the two connection 
10 pipes equal to the number of particles corresponding to 

a voltage of the voltage source. 

14. An electric network simulating apparatus 
according to one of claims 9 and 10, wherein the 
setting means includes 

15 means for, when a given one of the circuit 

elements has an impedance characteristic discontinu- 
ously changing, preparing a plurality of rules for the 
element cell for expressing the given circuit element 
and selecting one of the plurality of rules in 

2 0 accordance with the state of the connection pipe 

connected to the element cell. 

15. An electric network simulating apparatus 
according to claim 13, wherein the transferring means 
and the simulating means include 

25 means for simulating the state of each element 

cell at initial time so as to simulate a transient 
phenomenon of the given circuit element having 
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nonlinearity as a function of time, simulating a 
behavior of the nonlinear element at an operating point 
advancing by a shortest time interval, by changing each 
parameter of a combination of the element cells having 
5 functions equivalent to the element cells, and 

simulating the transient phenomenon by repeating the 
change in parameter every time the shortest time 
interval has elapsed. 

16. An electric network simulating apparatus 

10 according to claim 14, wherein the transferring means 

and the simulating means include 

means for simulating a behavior of each element 
cell at initial time so as to simulate a transient 
phenomenon of the given circuit element having the 

15 impedance characteristic discontinuous ly changing, 

simulating a behavior of the nonlinear element at 
an operating point advancing a shortest time interval 
by executing the transferring means in accordance with 
the rule selected in accordance with the state of the 

2 0 connection pipe connected to the element cell, and 

simulating the transient phenomenon by repeating the 
simulating means every time the shortest time interval 
has elapsed. 

17. A storage medium storing a simulation program 
2 5 loaded and activated in a computer device, the program 

activating the computer device to generate: 

means for defining element cells representing 
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electric functions of a plurality of circuit elements 
and connection pipes representing wiring lines for 
connecting the circuit elements, defining an electric 
network current as the number of particles moving in 
5 the connection pipe per unit time, and defining an 

electric network voltage as the number of particles 
present in the connection pipe; 

means for setting beforehand, on the basis of 
definitions in the defining step, in units of element 
10 cells, a rule for expressing an electric function of 

each of the circuit elements in accordance with a state 
of the connection pipe connected to each of the element 
cells ; 

means for transferring particles between the 
15 element cell and the connection pipe in accordance with 

the set rule; and 

means for simulating the state of the electric 
network by updating the number of particles passing 
through a given connection pipe per unit time in the 
2 0 transferring step and the number of particles present 

in the given connection pipe, and performing the 
transfer means and updating process at least once. 

18. A storage medium storing a simulation program 
loaded and activated in a computer device, the program 
2 5 activating the computer device to generate: 

means for, after setting element cells represent- 
ing electric functions of a plurality of circuit 
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elements, intersection cells representing functions of 
electric wiring intersections, and connection pipes 
representing connections between the element cells 
and the intersection cells, defining a current of 
5 an electric network as the number of particles moving 

in the connection pipe per unit time, and defining 
a voltage of the electric network as the number of 
particles present in the connection pipe; 

means for setting beforehand, on the basis of the 

10 definitions in the defining step, in units of element 

cells, a rule expressing an electric function of each 
of the circuit elements in accordance with a state of 
the connection pipe connected to the element cell, and 
setting beforehand, in units of intersection cells, 

15 a rule so that the numbers of particles present in the 

connection pipes connected to the intersection cell 
are equal to each other and a sum of the numbers of 
particles transferred at the intersection cell becomes 
zero; 

2 0 means for transferring particles between the 

element cell and the connection pipe and between the 
intersection cell and the connection pipe on the basis 
of the rules set in the setting step; and 

means for simulating the state of the electric 

25 network by updating the number of particles passing 

through a given connection pipe per unit time and the 
number of particles present in the given connection 
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pipe in the transferring means and performing transfer 
and updating processes at least once. 

19. A storage medium storing the simulation 
program according to one of claims 17 and 18, wherein 

5 the setting means includes 

means for setting, when a given one of the circuit 
elements is a current source, a rule for extracting the 
number of particles corresponding to a current value 
per unit time from one of two connection pipes 
10 connected to an element cell expressing the given 

circuit element and giving the number of particles 
equal in number to the number of extracted particles 
to the other one of the two connection pipes. 

20. A storage medium storing the simulation 

15 program according to one of claims 17 and 18, wherein 

the setting means includes 

means for setting, when a given one of the circuit 
elements is a voltage source, a rule for making a 
difference between the number of particles in one of 

2 0 two connection pipes connected to an element cell 

expressing the given circuit element and the number of 
particles in the other one of the two connection pipes 
equal to the number of particles corresponding to a 
voltage of the voltage source. 

2 5 21. A storage medium storing the simulation 

program according to one of claims 17 and 18, wherein 
the defining means includes 



means for defining, when a given one of the 
element cells has nonlinearity as a function of time, 
the given circuit element as a combination of an 
element cell for a resistive element and one of an 
element cell expressing a current source and an element 
cell expressing a voltage source, the combination 
expresses linearity equivalent to a behavior of the 
given circuit element at given time; and 

the setting means includes 

means for setting, when a certain one of the 
circuit elements is a current source, a rule for 
extracting the number of particles corresponding to a 
current value per unit time from one of two connection 
pipes connected to an element cell expressing the 
certain circuit element and giving the number of 
particles equal in number to the number of extracted 
particles to the other one of the two connection pipes, 
and 

means for setting, when a specific one of the 
circuit elements is a voltage source, a rule for making 
a difference between the number of particles in one 
of two connection pipes connected to an element cell 
expressing the specific circuit element and the number 
of particles in the other one of the two connection 
pipes equal to the number of particles corresponding to 
a voltage of the voltage source. 

22. A storage medium storing the simulation 
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program according to one of claims 17 and 18, wherein 
the setting means includes 

means for, when a given one of the circuit 
elements has an impedance characteristic 
5 discontinuously changing, preparing a plurality of 

rules for the element cell for expressing the given 
circuit element and selecting one of the plurality of 
rules in accordance with the state of the connection 
pipe connected to the element cell. 

10 2 3. A storage medium storing the simulation 

program according to claim 22, wherein the transferring 
means and the simulating means include 

means for simulating the state of each element 
cell at initial time so as to simulate a transient 

15 phenomenon of the given circuit element having nonlin- 

earity as a function of time, simulating a behavior of 
the nonlinear element at an operating point advancing 
by a shortest time interval, by changing each parameter 
of a combination of the element cells having functions 

2 0 equivalent to the element cells, and simulating the 

transient phenomenon by repeating the change in 
parameter every time the shortest time interval has 
elapsed . 

24. A storage medium storing the simulation 
25 program according to claim 23, wherein the transferring 

means and the simulating means include 

means for simulating a behavior of each element 
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cell at initial time so as to simulate a transient 
phenomenon of the given circuit element having the 
impedance characteristic discontinuously changing, 
simulating a behavior of the nonlinear element at 
5 an operating point advancing a shortest time interval 

by executing the transferring step in accordance with 
the rule selected in accordance with the state of the 
connection pipe connected to the element cell, and 
simulating the transient phenomenon by repeating the 
10 simulating steps every time the shortest time interval 

has elapsed. 
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ABSTRACT OF THE DISCLOSURE 
According to an electric network simulating 
method, element cells representing electric functions 
of a plurality of circuit elements and connection pipes 
5 representing wiring lines for connecting the circuit 

elements are defined. A current is defined as the 
number of particles moving through the connection pipe 
per unit time, and a voltage is defined as the number 
of particles present in the connection pipe. A rule 
10 for expressing the electric function of each circuit 

element in accordance the state of the connection 
pipe is set beforehand in units of element cells. 
The particles are transferred between the element cell 
and the connection pipe in accordance with the rule. 
15 The state of the electric network is simulated in 

accordance with the number of particles passing through 
the connection pipe per unit time and the number of 
particles present in the connection pipe. 
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Panod=0. 000019 



Ct = 
C2 = 
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C4 = 
C5 = 
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R2 = 
katu 



3. 8E-08 
0. 000046 
6. 8E-09 
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4.3E-09 
0.0012 
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: 20 

= 0.65 
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VI = 200 

L1 = 0.0012 
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OUTPUT SECTION 



V(11)_max : 
V(11)_min : 
V(30)_max = 
V(30)_min : 
I (R1)_max : 
I (R1)_min : 
I (R2)_max 
I (R2)_min 
I (L1)_max 
I (L1)_min 
I (L2)_max 
I (L2)_min 



683.3895 
-1.087639 
278. 2092 
-388.625 
0.4600302 
-0. 6466231 
0. 5038309 
—0.563192 
1 . 837644 
-1.170115 
0.6056873 
-0.9513868 
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CELL LIST 


CELL NUMBER 


CONNE 
NODE 


CTION 
NUMBER 


UNIQUE INFORMATION 


NME$(1)="V1": 

NME$(2)="D1": 

NME$(3)="X1": 

NME$(4)="X2": 

NME$(5)="Y1": 

NME$(6)="R C2 ": 

NME$(7)="VC2": 

NME$(8)="X1": 

NME$(9)="RC3": 
NME$(10)="V C3 ": 

NME$(11)="V2": 
NME$(12)="S1": 

NME$(13)="X3": 
NME$(14)="P1": 

NME$(15)="X4": 
NME$(16)="X5": 

NME$(17)="RTPI": 
NME$(18)="R]- M |": 

NME$(19)="ITPI": 
NME$(20)="X6": 
NME$(21)="X7": 
NME$(22)="B rS |": 

NME$(23)="RTMI" 

NME$(24)="I TS |": 

NME$(25)="Y2": 

NME$(26)="R1": 

NME$(27)="R2": 

NME$(28)="X8": 

NME$(29)="RC5": 
NME$(30)="V C5 ": 

NME$(31)="RC4": 
NME$(32)="V C4 ": 
new_i=32 


NP(1)=l: 

NP(2)=1: 

NP(3)=2: 
NP(4)=9: 
NP(5)=0: 
NP(6)=3: 
NP(7)=6: 

NP(8)=2: 

NP(9)=5: 

NP(10)=8 

NP(11)=13: 

NP(12)=11: 

NP(13)=13: 
NP(14)=39: 

NP(15)=4: 
NP(16)=12: 
NP(17)=18: 
NP(18)=23: 

NP(19)=20: 
NP(20)=23: 
NP(21)=26: 
NP(22)=24: 

NP(23)=19: 

NP(24)=28: 
NP(25)=30: 
NP(26)=32: 
NP(27)=31: 
NP (28) -35 " 
NP(29)=33: 
NP(30)=34: 
NP(31)=29: 
NP(32)=38: 


NM(1)=0: 

NM(2)=2: 

NM(3)=3: 

NM(4)=10: 

NM(5)=7: 

NM(6)=6: 

NM(7)=7: 

NM(8)=3: 

NM(9)=8: 

NM(1)=9: 

NM(11)=10: 

NM(12)=15: 

NM(13)=14: 
NM(14)=15: 

NM(15)=17: 
NM(16)=20: 
NM(17)=2l: 
NM(18)=26: 

NM(19)=17: 
NM(20)=24: 
NM(21)=27: 
NM(22)=27: 

NM(23)=22: 

NM(24)=25: 
NM(25)=31: 
NM(26)=36: 
NM(27)=33: 
NMf28)-3fi" 
NM(29)=34: 
NM(30)=35: 
NM(31)=38: 
NM(32)=37: 


DTA(1)=V1 

N3(3)=4: N4(3)=5 
N3(4)=11: N4(4)=12 
N3(5)=14: 
DTA(6)=dt/C2 
DTA(7)=0 

N3(8)=4: N4(8)=5 

DTA(9)=dt/C3 

DTA(10)=0 

N3(12)=39I N4(12)=15: 
DAT(12)=1 

N3(13)=15: N4(13)=16 
DAT(14)=50: DAT2(14)=0: 
DAT3(14)=223 
N3(15)=18: N4(15)=19 
N3(16)=21: N4(16)=22 
DAT(17)=L1 * (1-ketu*ketu)/dt 
DAT (18) = (LI *L2-M*M)/ 

(M*dt) 

DAT(19)=0 

N3(20)=25: N4(20)=30 
N3(21)=28: N4(21)=29 
DAT(22)=L2* (1-ketu* ketu)/ 
dt 

DAT(23) = (L1 *L2-M*M)/ 

(M*dt) 

DTA(24)=0 
N3(25)=32 
DAT(26)=R1 
DAT(27)=R2 

v*0\c.OJ—of. IM*4 V^iO^ = I D 

DTA(29)=dt/C5 
DTA(30)=0 
DTA(31)=dt/C4 
DTA(32)=0 
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